Spontaneous Emission from a Fractal Vacuum 
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Spontaneous emission of a quantum emitter coupled to a QED vacuum with a deterministic 
fractal structure of its spectrum is considered. We show that the decay probability does not follow 
a Wigner-Weisskopf exponential decrease but rather an overall power law behavior with a rich 
oscillatory structure, both depending on the local fractal properties of the vacuum spectrum. These 
results are obtained by giving first a general perturbative derivation for short times. Then we 
propose a simplified model which retains the main features of a fractal spectrum to establish analytic 
expressions valid for all time scales. Finally, we discuss the case of a Fibonacci cavity and its 
experimental relevance to observe these results. 
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Spontaneous emission results from the coupling of a 
quantum system with a discrete energy spectrum (an 
"atom") to a quantum vacuum. This is an important 
and widely studied phenomenon both from fundamen- 
tal and applied points of view [lj. Spontaneous emission 
allows to probe properties of the quantum vacuum, its 
dynamics and correlations. The wide zoology of behav- 
iors of spontaneous emission depends on spectral prop- 
erties of the vacuum and on its coupling to the atom 
. A standard textbook description 0, [J considers the 
coupling to a vacuum having a smooth and non-singular 
density of modes of photons, in which case, the prob- 
ability for spontaneous emission follows the well-known 
Wigner-Weisskopf decay law, 



\Ue(t)f 



-r e (w e )t 



(1) 



Relevant definitions of the quantum amplitude U e (t) and 
of the inverse lifetime T e (ui e ) will be given below. This 
description has been further developed towards quantum 
emitters coupled to more complicated environments such 
as semiconductors, QED-cavities, photonic crystals and 
micro-cavities The existence of singularities in the 
spectrum of the vacuum leads to a qualitatively different 
behavior which has been studied in various cases Q ■ 

In this letter, we address the problem of spontaneous 
emission from an atom coupled to a vacuum whose spec- 
trum is characterized by a discrete scaling symmetry ex- 
pressed by the property 



fx(u) + Auj) — fi(oj) 



At (T(w + Aw)) -fi (T (w)) 



(2) 



where fx (ui) is the integrated density of modes, or spectral 
measure, and the dimensionless scaling parameter a and 
the map T(lu) provide a full characterization of the spe- 
cific discrete scaling symmetry Introducing J\f Uu (ui) = 
fx(ui)—fx (oj u )i the scaling relation ([2]) can be written more 
concisely as A/L u (w) = \Nt{w v ) (T (w)). A spectrum de- 
scribed by ((2]) is often called fractal (see [5( and refs. 
therein), a denomination that we shall retain all over but 



which covers a broad class of systems that extends be- 
yond the conventional self-similar character associated to 
fractals, such as the popular Cantor set. More involved 
examples include cavities made out of quasi-periodic, e.g. 
Fibonacci, sequences of dielectric layers, or cavities gen- 
erated from deterministic self-similar fractal structures 
like a Sierpinski gasket. The prominent feature of frac- 
tal spectra is that they are highly lacunar, possessing an 
infinity of gaps appearing at all scales. This is a direct 
consequence of their discrete scaling symmetry (|2|). It 
is worth noticing that for a singular spectrum, the in- 
tegrated density of modes fj,(ui) is usually well defined, 
while its derivative, the density of modes p(ui) = , 
may not be. 

Singularities of the spectrum satisfying ([2]) correspond 
to fixed points of the map T (ui). Around a fixed point uj u 
the map can be linearized, T(ui) ~ u) u + T'(uj u )(u) — uj u ). 
Thus, the spectral measure obeys the equation M Uu (ui) = 
JV^,, (ui u + T" (u) u ) (u> — u> u )), whose general solution is 



6] 



Af Uu (u) = \u,-u, u \ a F 
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(3) 



wher 



In a 



is the local (w u -depcndent) spec- 



ln|T'(w„)| 

tral exponent and T(x) is a periodic function of period 
unity [7[. Generically, a changes between zero and unity 
spanning the range between smooth continuous and point 
spectrum 

We show below that the coupling of an atom whose res- 
onance frequency is close to ui u leads to a similar scaling 
behavior of the time dependent decay amplitude U e (t), 
given in the long time limit by 



\Ue(t)\ 



2 = t - 2l 



mi 



(4) 



where Q(x + 1) = Q(x) is another periodic function and 
7 is a function of the spectral exponent a Q . The expo- 
nent 7, the real parameter A = In \T'(u> u ) \ and the specific 
form of the function Q(x) are direct consequences of the 
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specific scaling relation ((3]). A similar scaling has been 
obtained for the thermodynamics of a fractal blackbody 
[Toj . Thus, for a fractal spectrum, spontaneous emission 
does not follow the Wigner-Weisskopf exponential decay 
(H|), a result which could be partly anticipated from the 
singular nature of the spectrum. It is indeed known that 
the existence of spectral singularities leads to an alge- 
braic time decrease of the decay probability qualitatively 
different from the exponential behavior or its small time 
approximation given by the Fermi golden rule Q. But 
the existence of log-periodic fluctuations described by the 
function Q (x) is a direct consequence of the discrete scal- 
ing symmetry ([2]). 

Expression (HJ constitutes the main result of this letter. 
To establish it, we shall first recall some basic definitions 
and results. Then, we will give a general perturbativc 
derivation starting from the Fermi golden rule, essen- 
tially limited to small times. To go beyond this limit, 
we will consider a model general enough to include all 
relevant characteristics of a fractal spectrum as given by 
([2]), yet simple enough to allow for a thorough analyt- 
ical derivation. Finally, we will discuss the case of the 
Fibonacci cavity which provides a relevant example of 
quantum vacuum with discrete scaling symmetry which 
can be realized experimentally. 

A two- level quantum system (\g), |e)) whose Hamilto- 
nian is H e = huj e |e) (e|, with respective energies and 
hu! e , is coupled to the EM field described by Hp = 
Wfeaj^a/j, where k stands for an appropriate set of 
quantum numbers. The atom-photon interaction is de- 
scribed within the rotating wave approximation by the 

Hamiltonian H int = J2k i^k a \ \d) ( e l + h.c^j where only 
resonant and near-resonant transitions are considered. 
The matrix element Vk denotes the strength of the cou- 
pling, in general dependent on the atom's position. In 
the initial state |e,0fc), the atom is in the excited state 
and no photon is present. We aim to evaluate the prob- 
ability amplitude U e (t) to find the quantum system in 
the initial state a time t after it evolves with the total 
Hamiltonian H = H e + Hf + ffmt. This amplitude is 
defined by 



U e {t) = (e,0 k \U(t,0)\e,0 k ) 



(5) 



The evolution operator U(t, 0), written in terms of the 
resolvent operator G(z) = l/(z — H), is Q 



U(t,0) 
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dEe- lEt ' h {G + {E) - GL (£)). (6) 



The matrix element G e (z) = (e, 0fc|G(z)|e, Ok) = \j{z — 
E e — S e (z)) of the resolvent defines the self-energy S e (w± 
i0 + ) = A e (ui) =F ^T e (uj) in terms of two spectral func- 
tions, A e (w) and T e (uj) which respectively account for 
the shift and the spectral width of the atomic energy E e . 
The spectral function T e (ui) is related to the vacuum 



response function 



(7) 



which can also be expressed using the spectral measure 
p(u>) as 



Within the 
£ k (r) (e z ■ £ fe ) d ge , 
correlation function 
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dipolc 
, the 
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approximation, 
response $ e (i) is 



(8) 

V k 

the time 



* e (t) = ft~ 2 \d 9 e\ 2 (0 fc |4 (r, t)E\ (r, 0) |0fc) 



(9) 



of the electric field component E z (r, t) along the polariza- 
tion direction z of the atom. Here, d ge and r are respec- 
tively the dipole matrix element and the position of the 
atom. A convenient form of the probability amplitude 
U e (t) is given in terms of the Laplace transform $ e (s) of 
the correlation function Q : 



U e (t) 



1 

2~7ri 



a+ioo 



ds- 



s + $ e (s - ioj e ) 



(10) 



Previous definitions and results summarize the state of 
the art and allow to identify two relevant energy scales 
for the problem of spontaneous emission. One, T e (uj e ) ~ 
|I4| 2 /9(w e ), is the strength of the coupling between the 
emitter and the vacuum as defined from the Hamiltonian 
Hi n t. The second energy scale A is given by the spec- 
tral width of r e (w). Relation between the two allows to 
identify a dimensionless coupling parameter 



r e (w e ) 



(ii) 



In the weak coupling limit, j« 1, the quantum ampli- 
tude U e (t) is determined by the pole in (fTU)) . given by 



the approximate solution s 



■i(-we) - -|r e ( 



This leads straightforwardly to the Wigner-Weisskopf ex- 
ponential decay (JXJ) of the probability amplitude. At ex- 
tremely long times, t ^> r~ 1 (oj e ), this pole approxima- 
tion is not valid anymore, even in free space, and both 
the probability amplitude U e (t) and the correlation func- 
tion ( i ) e (i), are dominated by the singularity at the edge 
lu = of the vacuum spectrum Q. For an atom cou- 
pled to the ci-dimensional scalar QED vacuum, T e (uj) ~ 
p(ui) \£ (w)| 2 , where the density of states p(uj) ~ a/" 1 
and the amplitude of the electric field £ (u>) ~ yfuj, so 
that T e ~ uj d and $ e (i) ~ l/t d+1 , and according to (fT0|) . 
U e (t) ~ l/t d+1 . Such singularities have also been studied 
in the lifetime of electronic states in quantum mesoscopic 
metals (llj . In more structured vacuum spectra such as 
photonic crystals, the spectral function T e (ui) exhibits 
singularities of the type T e (u>) = C\u> — u> u \ a ^ 1 8(uj — uj u ) 



3 



around certain frequencies co u In that case the cou- 
pling is strong and the pole approximation is not ap- 
plicable. Instead, one obtains a generalized exponential 
decay, limited to small times well accounted for by the 
Fermi golden rule, U e (t) — 1 ~ Ct 2 ~ a (see also discussion 
below). At large times, it turns into an algebraic decrease 
U e (t) ~ t~( 2 ~ a \ possibly coexisting with a non-decaying 
component [J]. 

We now consider the coupling to a fractal vacuum 
spectrum obeying the scaling in ((3|). In this case the 
weak coupling limit becomes ill-defined, since the spec- 
tral function T e (uj) ~ p(uj) is a singular function with 
a vanishing variation scale A. Thus, according to (fTT|l . 
we are effectively in a strong coupling regime j » 1, 
even for a finite and weak coupling Vk- On the other 
hand, the short time perturbative approach remains ap- 
plicable. We start with this approximation to present 
an intuitively clear explanation of the effect of the vac- 
uum spectrum fractality on the decay dynamics. At short 
times one obtains from (1101) 



We{t)\' 



1 - 



dt'r e (t r ), 



(12) 



where using (J7J and the explicit expression of the corre- 
lation function $ e (t) we have 



r.W - f 



dn(k)\V k \ 



2 sin (w fe 



LO e ) t 



(13) 



(co k -w e )t 

While this expression is textbook materials, we wish to 
re-examine it in the context of a fractal spectrum con- 
sidered here. The sine function on the r.h.s of (p~3|) indi- 
cates that r e (t) is not the Fourier transform of the spec- 
tral function r e (w), but it is rather a wavelet transform. 
Generally, the wavelet transform S w (a,b) of a function 
s(x) is defined by [l2[ 



S w (a,b) = — [ dxs(x)w 
a J 



(14) 



It can be viewed as a mathematical microscope which 
probes the function s(x) at a point b with a magnification 
1/a and an optics specified by the choice of the specific 
wavelet w(x). Thus, T e (t) in fp~3|) is the wavelet trans- 
form of the spectral measure /j.(u>) at a frequency uj e , with 
a magnification t and w(x) = sinc(x) as a probe. An im- 
portant property of the wavelet transform is that it pre- 
serves a discrete scaling symmetry (|2} of the probed func- 
tion, here the spectral measure p{w) weighted by \Vk\ ■ 
For a smooth and continuous spectrum, the sine function 
probes the spectrum on energy scales of the order of t -1 
and in the long time limit it goes over 5(u>k ~ 0J e ) so that 
r e (i) becomes t- independent as expected. However there 
is no such well-defined limit for a fractal spectrum and 
inserting J3J into (I13p for cu e = ui u , we obtain instead [T3| 



T e (t) = t 1 - a f 1 



lnt 



In|T'(w„)| 



(15) 




FIG. 1: (Color online) Numerical results for T e (t), given by 
(|13p . for Fibonacci dielectric cavities Sj [defined in text] of dif- 
ferent lengths N. The photonic spectrum has been obtained 
from a simplified diagonal tight-binding (TB) description of 
the Fibonacci potential, which however exhibits the relevant 
fractal properties of a dielectric superlattice. The matrix ele- 
ment |F fc | 2 in was taken constant, which is equivalent to 
averaging T e (t) over the atom's position in the cavity. The 
dimensionless time t is in units of the TB hopping constant. 



and 



1 -f 



2-a 



lnt 



ln|T'(w u )| 



(16) 



where T\^(x) are periodic functions of their argument of 
period unity. This constitutes a short-time counterpart 
to the asymptotic result The behavior of r e (t) is 
illustrated in Fig. Q] for the specific case of a Fibonacci 
quasi-periodic dielectric cavity (to be discussed later on) , 
whose spectral measure admits locally a discrete scal- 
ing symmetry (j3|). We observe, as predicted by (fl~5j) . an 
overall power law behavior with a ~ 0.8, explained by 
the renormalization group analysis of [14|, and also log- 
periodic oscillations around it, which are the fingerprint 
of the underlying fractal structure of the spectrum. Note 
that the log-periodic oscillations characteristic of a frac- 
tal spectrum are already noticeable for systems of finite 
size and remain unchanged as length increases. 

The previous derivation is based on the Fermi golden 
rule valid at short times. In order to go beyond this 
regime, a more elaborate derivation is needed. To that 
purpose, we consider the following model for the spectral 
function: 



r e (w) = 



c 



ir \uj - U! u \ 



1 + A cos 



2tt 

T 



In- 



, (17) 

where < a < 1 is the local spectral exponent intro- 
duced in ((U, C is the coupling strength, and < A < 1 
and n define the modulation amplitude and phase re- 
spectively. For simplicity, we let — oo < ui < oo. This ex- 
pression exhibits the basic features required to describe 
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a fractal structure of the vacuum as denned in ([2J and 
(J3]) with A = In \T' (uj u )\. We have approximated the log- 
periodic function by its first harmonic, which happens to 
be a good approximation as shown in related situations 
ljj 16 1 . In the absence of log-periodic modulation, i.e. 
for A — 0, we recover the case of a singularity in the spec- 
trum similar to that already studied in the literature Q ■ 
The important point here is that the modulation results 
from the specific scaling properties of the spectrum de- 
fined in ((21) and discussed subsequently. The model (fl7|) 
allows to obtain closed analytical expressions of the quan- 
tum amplitude U e (t) and therefore to study its behavior 
both in the short time limit (in a way very similar to the 
discussion which follows (|12j) ) and in the long time limit 
which was not possible starting from the Fermi golden 
rule. From ([7]), we obtain 



*e(t) = C 



2e 



irt a 



F(a) cos 



ira 



AlmF(t) 



(18) 



where T(x) is the Euler Gamma-function and we have de- 
fined F(t) = {nt) 2l7r/x cosh (2f + i^s) r (a - 2jE). We 

thus immediately notice that <& e (i) is not short range in 
time. Consequently, there is no Wigner-Weisskopf expo- 
nential decay as in ([lj. While this is true already for 
A = 0, note that the log-periodic modulation of r e (w) 
adds to $ e (i) an oscillatory log-periodic term, which fur- 
ther modifies the behavior of U e (t). To investigate this 
point in more detail, we need to study the pole structure 
in (fTO)) . The Laplace transform of $ e (i) is 



C 



csc ™ + — [F(z) -F*(z 
2 2i X y ' V ' 



where z 



(19) 

= s + iAui, Auj = u) e — u) u and F(z) = 
-> sinlT 1 - i^j. The poles s n in (JTUJ are 

solutions of s + $ e (z) = 0. For the sake of simplicity, 
we consider Au; = 0, in which case all the poles come in 
complex conjugate pairs. For A not too small, namely 
A > esc (tP) / A , the poles are infinitely dense near 
the origin and are distributed log-periodically with the 
distance to it. In particular, for |s n | 2 ~ a <C C, their ex- 
pression in the lower half-plane is [171 ] 



(20) 



with 8o = y~ In (\A\ sin Tp) and n is an integer. The 
corresponding residues acquire a rather simple form, 
res(s n ) rj — sin (^) 2^js^ Q e s "*. Then, assuming 
e -27r 2 /A ^ i ^ ne ]-, ranc h cut contribution to (TTU1) can 



be neglected, and the probability amplitude U e (t) 



e iUet U e (t), in the long time limit C(|sin6» |0 2 " > 1 



can be written as 1181 



U e (t) 



Asin^p 



+ 00 

E • 

n— — oo 



(21) 
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FIG. 2: (Color online) Amplitude U e (t) as a function of time 
for the model ([T7J) with a = |, A = 1 and A = (blue 

thick line) or A — 1 (red thin line). Inset shows f7 e (£) in a 
semi-log scale. 



Here, the upper summation limit was extended to +oo 
exploiting the condition on t. It is possible to make 
further approximations and show that only a few terms 
in the sum give the main contribution at a given time. 
Therefore, locally, the decay envelope appears as a slow 
exponential, while it is algebraic at large time scales, 

A sin(-7ra/2) 



Pe(t)\ 



C(\sin9 \ty 



discrete scaling symmetry displayed in (|2Tj) . namely 
U e {pt) = P a ~ 2 U e (t) with/3 = e\ 



17]. For comparison, the decay 

asymptotics for A = is |E/e(t)| — 1/Ct 2 ~ a , as deter- 
mined by the branch cut in (|T0| . More interesting is the 



22) 



which is straightforward when using (|20[) in (|2Tj). More 
generally, if for large t one can neglect the free term s in 
the denominator in (fTtJl) . then ([2"2"]) immediately follows 
from scaling properties of <fr e (s — ioj e ). Thus, as a result 
of ([2|), U e (t) takes precisely the form ((4|) with 7 = 2 — 
a. The log-periodic function Q can be calculated from 
(|2T|). The decay amplitude U e (t) is plotted in Fig. H 
This figure makes explicit that U e (t) has a much slower 
decay for A = 1, than for A = 0, i.e., without the log- 
periodic modulation. As noted above, the decay envelope 
appears as exponential since the power law in |4]) shows 
up at larger time scales than displayed in Fig. [2l The 
oscillating structure of U e (t) shows beats, which result 
from interferences between the few dominant terms in 

A physical realization of the previous results is pro- 
vided by an atom, or a quantum dot, coupled to a vac- 
uum field of a quasi-periodic, e.g. Fibonacci, dielectric 
cavity [H, [20| . The latter is made of a sequence of 
slabs of two types, A and B, of width and refractive 
index denoted respectively by <Ia , cLb and tla , jib . Fi- 
bonacci sequences Sj of such slabs are constructed by 
recursion Sj>3 = [Sj-2Sj-i], with Si = B, S% = A. 
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The spectral properties and the spatial behavior of the 
corresponding eigenmodes in Fibonacci and some other 
quasi-periodic structures have been extensively studied 
2l|. The spectrum of the Fibonacci system Sj is highly 
fragmented, with the degree of fragmentation increasing 



with j and the contrast tla/tib [IJ, [22|. More specif- 
ically, the spectral measure h(uj) has a discrete scaling 
symmetry ([2]) governed by specific p-cycles of the as- 
sociated renormalization group transformation, whereas 
a = a v with a = ■ Near fixed points of a given 

p-cycle map T (uj) can be linearized, and /x (w) obeys 
® with a = ^ = ^f- and A = ln\T' For 
instance, for d^riA — dsns = d, fixed points of a 6- 
cycle are determined by = 2 " 4 t " 1 , where c is the 
vacuum speed of light and n = 0,1,.... Linearizing 
the map at these points, one obtains e A = |T'(w n )| = 

1 + 8ry 4 + 4t)Vl + ¥ with r? = A + 2a) 0. Ex- 
pcrimcntally, a reasonable contrast can be achieved with 
n A = 1.45. and n B = 2.23 [H, so that a = 0.898 and 
A = 3.2 



23|. 



In conclusion, we have argued that spontaneous emis- 
sion of a quantum emitter coupled to a vacuum char- 
acterized by a discrete scaling symmetry exhibits un- 
usual and measurable features. We have shown that the 
quantum probability amplitude never follows the Wigner- 
Weisskopf decay law (JlJ, but rather an overall power 
law whose exponent depends on the fractal properties 
of the vacuum. Moreover, we have shown that the de- 
cay probability is characterized by log-periodic fluctua- 
tions, both at short and long times. These fluctuations 
constitute an unambiguous fingerprint of the underlying 
fractal structure. Underneath the slow algebraic decay 
and log-periodic modulation, the probability amplitude 
was shown to exhibit a rich dynamics with nearly peri- 
odic oscillations and beats resulting from the interference 
between the main contributions to the pole structure. 
Finally, we have discussed the example of a Fibonacci 
cavity, a non fractal device but whose spectrum is char- 
acterized by a scaling symmetry ([2J , and we have argued 
that it may be a possible candidate to observe the un- 
usual features of spontaneous emission described in this 
letter. Spontaneous emission constitutes a specific way 
to probe a fractal nature of spectrum. Thermodynamics 
[l6( and spatial correlations [10( provide other examples 
of physical probes. Let us mention finally, that the ap- 
proach developed in this letter could be extended to re- 
lated problems in quantum mesoscopic physics (Coulomb 
blockade) and in the study of quantum vacuum effects 
such as the dynamical Casimir effect and the Schwinger 
and Unruh effects. 
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